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Abstract 

B functions are a class of relatively complicated exponentially decaying basis functions. Since the molecular 
multicenter integrals of the much simpler Slater- type functions are notoriously difficult, it is not at all 
obvious why B functions should offer any advantages. However, B functions have Fourier transforms 
of exceptional simplicity, which greatly simplifies many of their molecular multicenter integrals. This 
article discusses the historical development of B functions from the perspective of the interaction between 
mathematics and theoretical chemistry, which traditionally has not been very good. Nevertheless, future 
progress in theoretical chemistry depends very much on a fertile interaction with neighboring disciplines. 

Keywords: Electronic structure theory, exponentially decaying basis functions, B functions, multicenter 
integrals, interdisciplinary collaboration. 
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1 Introduction 

Molecular electronic structure theory tries to explain 
the nature of molecules and their properties by doing 
calculations on the basis of quantum mechanics. This 
approach is obviously possible, but as we all know, 
numerous obstacles have to be overcome. Already in 
1929, Dirac had stated [U p. 714]: 

The fundamental laws necessary for the 
mathematical treatment of a large part 
of physics and the whole of chemistry are 
thus completely known, and the difhculty 
lies only in the fact that application of 
these laws leads to equations that are too 
complex to be soluble. It therefore be- 
comes desirable that approximate prac- 
tical methods of applying quantum me- 
chanics should be developed, which can 
lead to an explanation of the main fea- 
tures of complex atomic systems without 
too much computation. 

It is noteworthy that Dirac only mentioned atomic 
calculations. Apparently it was at that time more or 
less unthinkable to do meaningful molecular calcu- 
lations. But since the time of Dirac, there has been 
enormous progress in particular in numerical and ap- 
plied mathematics. Moreover, we have witnessed a 
computer revolution which now gives us a comput- 
ing power that was beyond imagination in 1929. Ac- 
cordingly, we now have much better chances of ac- 
complishing our aims. Chemical insight via quantum 
mechanical calculations is no longer utopia. 

Nevertheless, the central methodological problem 
addressed by Dirac - the development of powerful 
approximation techniques - has not really changed. 
The reduction of chemical phenomena to manageable 
computational problems has been and still is far from 
trivial. In spite of the theoretical and technologi- 
cal advances of recent years, the computational com- 
plexity of sufficiently accurate quantum mechanical 
calculations is still prohibitive, and molecular elec- 
tronic structure calculations are only feasible if we 
make some partly very drastic approximations. 

Thus, progress depends quite crucially on our 
ability of reformulating and translating chemical in- 
tuition to a mathematical formalism, which does bet- 
ter than the currently available approximation pro- 
cedures and which ultimately helps to make more 
accurate calculations on larger molecules possible. 

Optimists might argue that future advances in 
computer hard- and software will be enough to 



achieve this goal. However, I am not so optimistic. 
In the long run, we will also need substantial progress 
in mathematics. But this will to a large extend de- 
pend on a fruitful collaboration between theoretical 
chemists interested in methodological questions and 
mathematicians interested in scientific applications. 

In theoretical chemistry, there is the widespread 
believe that the parentage of atoms facilitates our at- 
tempts of understanding the electronic structure of 
molecules. Thus, in molecular calculations on the ba- 
sis of the Hartree-Fock-Roothaan equations [2 [3l |4j , 
effective one-particle wave functions (MOs) are usu- 
ally approximated by a linear combination of so- 
called basis functions or atomic orbitals (AOs) cen- 
tered at the different nuclei of the molecule: 

n N 

V'W = E E C:^n{r-Ra). (1.1) 
Q— 1 n— 1 

The fundamental entities in this LCAO-MO ap- 
proach are the basis functions {'fin{r)}n=i- Basis 
functions and their derivatives up to a certain order 
must be square integrable with respect to an inte- 
gration over the whole R'^, and they also have to be 
complete in the corresponding Hilbert and Sobolev 
spaces, but otherwise they are essentially arbitrary. 

There are, however, two additional requirements 
of utmost practical relevance which good basis func- 
tions should satisfy: 

1. Already short linear combinations of the type 
of (jl.ip should provide sufficiently accurate ap- 
proximations, i.e., the summation limit N in 
(|1.1|) should be as small as possible. 

2. The molecular multicenter integrals, which oc- 
cur inevitably in the LCAO-MO approach, 
should possess a manageable complexity. 

It is not too difficult to satisfy the first require- 
ment by taking into account what we know about 
the exact solutions of effective one-particle atomic 
and molecular Schrodinger equations. Kato [5] had 
shown that the singularities of the potential of atomic 
and molecular Hamiltonians produce discontinuities 
of the wave functions commonly called Coulomb or 
correlation cusps. Moreover, exact solutions of these 
Hamiltonians decay exponentially at large distances 
(see for example [6j[7] and references therein). 

These properties of exact solutions clearly favor 
Slater- type functions [51 [S] , which were originally in- 
troduced to provide computationally convenient ap- 
proximations to numerically determined solutions of 
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effective one-particle equations. In unnormalizcd 
form, they can be expressed as follows: 



Shavitt [El Eq. (55) on p. 15] - defined by Filter and 
Steinborn fW\, Eq. (2.14)] as follows: 



x'jAf3,r) = (/3r)^-^-ie"'''-yf 



yl'iPr). (1.2) 



Here, r £ 



solid harmonic and Y^'^ 



{l3r)^Yl^{e,(j)) is a regular 
', (f) is a (surface) spherical 
harmonic, /3 > is a scaling parameter, is a gen- 
eralized principal quantum number which is in most 
cases, but not always a positive integer > L -I- 1, and 
L and M are the usual (orbital) angular momentum 
quantum numbers. 

Slater-type functions were used with considerable 
success in atomic calculation, and there is no rea- 
sonable doubt that their use in molecular calcula- 
tions would also be highly desirable. Unfortunately, 
Slater-type functions are not able to satisfy the sec- 
ond requirement mentioned above: The evaluation 
of their molecular multicenter integrals turned out 
to be extremely difficult, and in spite of the heroic 
efforts of numerous researchers, no completely satis- 
factory solution has been found yet. Actually, the ef- 
ficient and reliable evaluation of multicenter integrals 
of Slater-type functions is among the oldest math- 
ematical and computational problems of molecular 
electronic structure theory. A review of the older 
literature can be found in an article by Dalgarno 
[lOJ. Another useful review is the one by Harris and 
Michels JLlJ. 

Because of the problems with the efficient and 
reliable evaluation of the notorious multicenter in- 
tegrals of Slater- type functions, molecular electronic 
structure calculations are now dominated by Gaus- 
sian basis functions which have many disadvantages: 
They are unphysical because they cannot reproduce 
a Coulomb cusp, and for large distances they also 
decay much faster than exponentially. Accordingly, 
fairly large values of N in (|l.ip are needed to pro- 
duce sufficiently accurate results. In fact, Gaussian 
basis functions have only one, albeit decisive advan- 
tage: Their molecular multicenter integrals can be 
computed comparatively easily. 

In spite of the clear dominance of Gaussian 
functions, the search for alternative basis functions, 
that are more physical than Gaussians and whose 
multicenter integrals can be computed more easily 
than the notoriously difficult multicenter integrals of 
Slater-type functions, has been continued by numer- 
ous researchers. In this article, I want to discuss the 
history and most interesting features of a class of 
exponentially decaying basis functions, the so-called 
B functions, that were ~ based on previous work by 



S",(/3,r) - 



2"+^(n + £)! 



(1.3) 



Here, /3 > 0, n e Z, and A:„_i/2 is a reduced Bessel 
function. If K^{z) is a modified Bessel function of 
the second kind [Ml p. 66], the reduced Bessel func- 
tion with in general complex v and z is defined as 
follows [H Eqs. (3.1) and (3.2)] 



K{z) = (2A)i/2z''X,(z). 



(1.4) 



If the order v is half-integral, v = n + 1/2 with 
n G No, a substantial simplification takes place: 
Then, the reduced Bessel function is an exponential 
multiplied by a terminating confluent hypergeomet- 
ric series iFi (see for example [TBI Eq. (3.7)]): 

= 2"(l/2)„e-^f^i(-n;-2n;2z). (1.5) 

Obviously, B functions are fairly complicated mathe- 
matical objects. Moreover, (|1.3p and (|1.5p imply that 
B functions can be expressed as a linear combination 
of Slater- type functions defined by (|1.2p . Hence, it 
is not at all clear why the use of the comparatively 
complicated B functions should offer any advantages 
over Slater-type functions which possess an excep- 
tionally simple explicit expression. 

If we form finite linear combinations of Slater- 
type functions and do some mathematical manipu- 
lations, then the complexity of the resulting expres- 
sion normally increases, depending on the number of 
Slater-type functions occurring in the linear combi- 
nation. In fortunate cases, however, it may happen 
that most terms of the resulting expression cancel 
exactly. Thus, a significant reduction of complexity 
by forming linear combinations is also possible. 

The topic of this article arc B functions and their 
role in the theory of multicenter integrals. However, 
research on multicenter integrals is essentially math- 
ematical in nature. Quantum mechanics only deter- 
mines which integrals we evaluate, but the techniques 
employed for their evaluations are entirely mathe- 
matical. Physical insight is of secondary importance. 
What really matters are mathematical skills. 

Thus, a discussion of B functions immediately 
raises another, very profound question about the role 
of mathematics in molecular electronic structure the- 
ory, and how well and how rapidly new advances in 
mathematics are incorporated into electronic struc- 
ture theory. 



Ernst Joachim Weniger: The Strange History of B Functions . . . 



3: Shavitt's Gauss Integral Representation 



4 



Essentially this boils down to the question 
whether and how well theoretical chemists interested 
in methodological questions and mathematicians are 
interacting and - in the best of all worlds - even col- 
laborating. As I had argued in [TTl Section 1], science 
is now highly fragmented and consists of almost com- 
pletely disjoint subfields. As a consequences of this 
probably unavoidable specialization, communication 
between researchers belonging to different subfields 
has deteriorated quite a bit and no longer functions 
as well as it should. 

Molecular electronic structure theory is a highly 
interdisciplinary research topic whose progress de- 
pends crucially on a functioning communication with 
and input from neighboring disciplines. Ultimately, 
computational chemistry has become feasible be- 
cause of the spectacular advances in computer hard- 
and software. However, mathematics is the basis of 
all computing and not all problems can be solved 
by using more powerful computers. One should also 
not forget that Pople, whose suite of Gaussian pro- 
grams has contributed greatly to the computeriza- 
tion of chemistry, had done his undergraduate stud- 
ies in mathematics, and later he had been reader in 
mathematics at the University of Cambridge. 

It is my central hypothesis that it should be in the 
self-interest of theoretical chemists to improve com- 
munication with mathematicians. Traditionally, this 
communication has been fairly bad. In my opinion, 
the history of B functions provides ample evidence 
supporting this claim. 



2 Bessel Polynomials 

Polynomials associated with the names of Legendre, 
Gegenbauer, Jacobi, Hermite, and Laguerre certainly 
belong to the pillars of classical analysis, and they 
are also of considerable importance in quantum the- 
ory and other scientific disciplines. 

There are several other, less important classes of 
polynomials with well established properties. Many 
of these polynomials are completely ignored by the- 
oretical physicists or chemists, because they have no 
obvious physical relevance. Such a pragmatic atti- 
tude is certainly justified. Nevertheless, occasional 
surprises cannot be ruled out. 

For example, in the later stages of my PhD thesis 
[18] I became aware of a book by Grosswald [19] from 
which I learned that the polynomial part in (jl.Sp and 
thus also reduced Bessel functions had already been 
treated in the mathematical literature. For the poly- 



nomial part of reduced Bessel functions, the notation 



-1-1/2 



(^) 



(2.1) 



is used in the mathematical literature [191 Eq. (1) 
on p. 34]. Together with some other, closely related 
polynomials, the 9n{z) are called Bessel polynomials 
[19]. According to Grosswald [HI Section 14], they 
have been applied in such diverse and seemingly un- 
related fields as number theory, statistics, and the 
analysis of complex electrical networks. 

The Bessel polynomials 6'„(z) occur also in a com- 
pletely different mathematical context which later 
became very important for mc. In the book by Baker 
and Graves-Morris [20, p. 8] on Pade approximants, 
it is remarked that Pade had shown in his seminal 
thesis 2 1 that the Pade approximant [n/m] with 
n, m G No to the exponential function exp(2;) can be 
expressed as the ratio of two terminating confluent 
hypergeometric series [ID] Eq. (2.12)]: 



[n/m] 



iFi(—n; —n — m; z) 
L-Fi(— m; —n — m; —z) 



(2.2) 



Accordingly, the diagonal Pade approximant with 
n — m to the exponential function can be expressed 
as the ratio of two Bessel polynomials: 



[n/n] 



0n{z/2) 
0n{-z/2) ' 



n e No . 



(2.3) 



3 Shavitt's Gauss 
Representation 



Integral 



We can safely assume that Bessel polynomials were 
completely unknown among theoretical chemists. 
Thus, the usefulness of reduced Bessel functions had 
to be discovered differently. 

It seems that Shavitt was the first theoretical 
chemist who noticed that reduced Bessel functions 
have in spite of their undeniable complexity some 
very useful features. In tables of Laplace transforms, 
one can find the following relationship (see for exam- 
ple fIT, Eq. (13.42) on p. 338]): 

s^/^ K.{as^/^) = -J^^yJT^ ■ (3-1) 

Setting s = z"^ and a — \ yields the following remark- 
ably simple Gaussian integral representation of the 
reduced Bessel function [12, Eq. (55) on p. 15]: 

Kiz) ^ iV^y^' I dt. (3.2) 
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Shavitt defined the reduced Bessel function accord- 
ing to kx{z) = z^K\{z). Unfortunately, this nota- 
tion is misleading because k\{z) can easily be con- 
founded with the spherical modified Bessel function 
k^{z) = ['k/{2z)Y/'^K^+i/2{z) with n e No. There- 
fore, Steinborn and Filter ^ Eqs. (3.1) and (3.2)] 
defined reduced Bessel functions via (|1.4p . The in- 
troduction of the prefactor (2/7r)^/^ in (|1.4p has the 
additional advantage that it simplifies most formulas. 

A Gauss integral representation for Slater-type 
functions can also be constructed. The starting point 
would be the following following Laplace transform 
(see for example Eq. (5.94) on p. 259]): 



1/2^ = [2/7r]i/2 

iV(8t) 

-^D2..+i{a/[2tf'^)dt. (3.3) 



s" exp(^— as"' "j = [2/7r] 

foo g-st-a^/(8t) 



{2tY 



Here, D21/+1 (a/[2t]^/'^) is a parabolic cylinder func- 
tion [14, p. 324] which is a special case of the Whit- 
taker function W^,^{z) of the second kind [111 p. 
296]. Accordingly, the Laplace transform (|3.3p is 
much more complicated than (|3.ip . Obviously, this 
conclusion applies also to the Gauss transform for 
Slater- type fmictions that can be derived from (j3.3|) . 

The remarkably compact Gauss integral represen- 
tation ()3.2|) inspired Shavitt [121 P- 16] to propose the 
use of reduced Bessel functions as alternative expo- 
nentially decaying basis functions in electronic struc- 
ture calculations. This was done in articles by Bishop 
and Somorjai [23], Somorjai and Yue 24], and Yue 
and Somorjai [25j . 

The Gaussian function in (|3.2p makes it possi- 
ble to reduce multicenter integrals involving reduced 
Bessel functions to the corresponding integrals of 
Gaussians, which can be evaluated comparatively 
easily, multiplied by integrals over nonphysical vari- 
ables. Unfortunately, the remaining integrations can 
normally only be done numerically. Therefore, this 
approach is apparently not efficient enough to pro- 
vide a viable approach for the evaluation of multi- 
center integrals of reduced Bessel functions. 

4 The Gegenbauer Addition 
Theorem of Reduced Bessel 
Functions 

Multicenter integrals are difficult to evaluate because 
the integration variables are not separated. Princi- 
pal tools, which can accomplish such a separation 



of variables, are so-called addition theorems. These 
are expansions of a given function /(r ± r') with 
r,r' G in terms of other functions that only de- 
pend on either r or r'. The best known example 
of such an addition theorem is the Laplace expan- 
sion of the Coulomb potential in terms of spherical 
harmonics, which is nothing but the well known gen- 
crating function [l — 2xz + z"^] of the Legendre 
polynomials [TH p. 232] is disguise. 

The modified Bessel function w^'^ K,y{'-fw) with 

1 /2 

w — [p^ -t- — 2rpcos6'] , < p < r, and 
e C \ No possesses the following Gegenbauer-type 
addition theorem [H pp. 106 - 107]: 

w-^K^i-fw) = 2'^7~T(i/)(rp)-'' 

OQ 

X ^ C;:(cos0) /.+„(7P) K,,+nhr) . (4.1) 

n=0 

Here, C^(cos6') is a Gegenbauer polynomial [lH p. 
218], and Iy+n{ip) is a modified Bessel function of 
the first kind [3 p. 66]. 

Comparison with (|1.4p shows that the addition 
theorem (I4.ip is essentially a Gegenbauer-type addi- 
tion theorem of the reduced Bessel function k-y(^w). 
A very important special case occurs for = 1/2 
because then the Gegenbauer polynomial becomes a 

1 /2 

Legendre polynomial according to C„ (x) ~ Pn(x) 
[l4l p. 219], yielding the Legendre- type addition the- 
orem of the Yukawa potential [HI p. 218]: 



W 



{rp)-^'^ }_^(2n+l)P„(cos0) 

n=0 



^ In+l/2{lP)Kr,+i/2{ir) . (4.2) 

With the help of the so-called addition theorem of 
the Legendre polynomials [2S1 p. 303] 

47r 



Pf [cos d] = 



E 



2^+1 



[Yr(r/r)YYr{r'/r'), (4.3) 



where cos9 — rr'/{rr'), the Legendre-type addition 
theorem ()4.2p can be converted to an expansion in 
terms of spherical harmonics. 

This addition theorem for the Yukawa potential 
was the starting point for the so-called zeta function 
method of Barnett and Coulson [27 . They tried to 
construct an addition theorem for Slater-type func- 
tions by applying suitable generating differential op- 
erators to this addition theorem. Unfortunately, this 
idea did not lead to a complete success in the case 
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of Slater- type functions. Due to technical problems 
it was not possible to perform all differentiations in 
closed form. Consequently, the coefficients of the 
zeta function expansion could only be computed re- 
cursively, but not in explicit form [271 HH] ■ 

Later, numerous other researchers worked on ad- 
dition theorems for Slater-type functions (see for ex- 
ample Refs. 33 - 50, 52, 65, 84, 85, 90, 91]). But 
all these addition theorems turned out to be fairly 
complicated. I suspect that there is still plenty of 
room for improvements. 

The technical problems with the zcta function 
method of Barnett and Coulson [57| in special and 
with addition theorems for Slater-type functions in 
general inspired Steinborn and Filter to derive an 
addition theorem for reduced Bessel functions di- 
rectly from the Gegenbauer expansion (|4.ip . This 
was much easier than the derivation of analogous ad- 
dition theorems for scalar Slater-type functions. 

Let us assume that a Gegenbauer-type addition 
of the type of (14. ip for some function is known. If 
the Gegenbauer polynomials in this expansion can 
be replaced by a finite sum of Legendre polynomials, 
a rearrangement of the original expansion yields an 
expansion in terms of Legendre polynomials. Then, 
we only need l|4.3p to obtain an expansion in terms 
of spherical harmonics. 

The practical realization of this obvious idea was 
apparently not so easy. As discussed by Steinborn 
and Filter [15, pp. 269 - 270], many authors had 
quite a few problems with the determination of ex- 
plicit expressions for the coefficients of the expan- 
sion of Gegenbauer polynomials in terms of Legendre 
polynomials. Also Steinborn and Filter constructed 
very messy expressions for these coefficients which 
are restricted to certain superscripts of the Gegen- 
bauer polynomial [301 Section 3]. 

However, already at that time a much more con- 
venient expression for these expansion coefficients 
been available in the mathematical literature. In Ex- 
ercise 4 on p. 284 of Rainville's book [31], one finds 
the following expansion of Gegenbauer polynomials 
in terms of Legendre polynomials, but no explicit 
proof (compare also 02^ Eq. (5.2)]): 



Lm/2J 
s=0 

X (2m — 4s 



ifJ-)m-s ifJ- - 1/2)^ 
(3/2)„_,s! 



(4.4) 

Here, [m/2\ denotes the integral part of m/2. 

Expansion (j4.4p can be proved via the explicit ex- 
pression [33l Eq. 7.313.7 on p. 836] for the integral 



J^^{l-x)°'{l~zY-^/'^C^^{x)C^{x)dx. One only has 
to set v = 1/2 and perform the limit a 0, which 
requires, however, some algebraic trickery. 

In this way, the following addition theorem for a 
reduced Bessel function with half-integral order can 
be derived [32l Eq. (5.5)]): 



kn-i/2{P\r< ±r> 



(-l)"87r 



n-l/2fr, \n-l/2 



(/3r<)"-^/^(/3r>) 



(2n- 1)!! 

oo i 

E E iTlY[Yr{rJr^)YYr{ry/r^) 



=0 m=-e 



■J^ (-n)^(l/2 - n)i+^ 
^!(3/2),+. 



i£ + 2v-n+l/2) 



X L 



+ 1/2 (/??■< )-f^«+2i/-n+l/2 (/3r> ) . (4.5) 



£+2i/-n 

Here, we have |r<| < |r>|. The original version of 
this addition theorem as derived by Steinborn and 
Filter [151 Eq. (3.4)] still contains unspecified expan- 
sion coefficients of Gegenbauer polynomials in terms 
of Legendre polynomials. 

The addition theorem (|4.5p was quite consequen- 
tial for my later scientific interests. In my diploma 
thesis (34j . which was published in condensed form 
in [35j . I used this addition theorem for the evalua- 
tion of simple multicentcr integrals of reduced Bessel 
functions. 

5 Convolution Theorems of B 
Functions 

The original topic of Filter's PhD thesis [36] was the 
construction of addition theorems and their use for 
the evaluation of multicentcr integrals. During these 
investigations, it became obvious that the use of ad- 
dition theorems alone would not suffice to permit an 
efficient and reliable evaluation of the complicated 
multicentcr integrals of exponentially decaying func- 
tions. Thus, alternative approaches were necessary. 

Starting from the well known addition theorem 
of the so-called modified Helmholtz harmonics. Fil- 
ter succeeded in his PhD thesis [33 Sections 7 - 9] to 
derive with the help of some very skillful mathemat- 
ical manipulations remarkably compact expressions 
for convolution, nuclear attraction, and Coulomb in- 
tegrals of B functions which were later published in 
03. ,37,. 

Probably the most spectacular result is the ex- 
tremely compact expression for the overlap integral 
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of two B functions with equal scaling parameters 
P > 0, which is simply a finite linear combination 
oiB functions [13l Eq. (4.3)]: 



— ^ (4m2|^imi|£m2 - mi) 

A£ 



E<-.'(r) 

t=o ^ ^ 



-t+ 



,AP,R). (5.1) 



Here, (^2W2|£iTOi |^?7i2 — mi) is a so-called Gaunt 
coefficient, which is the integral of three spherical 
harmonics over the surface of the unit sphere in M'^. 
A compact review of the properties of Gaunt coef- 
ficients plus additional references can be found in 
[T71 Appendix C]. Moreover, it follows from the se- 
lection rules satisfied by the Gaunt coefficient (see 
for example [351 Eq. (3.1)]) that the summation limit 
A£ = (£i -I- £2 — £)/2 is a non-negative integer. 

For overlap integrals with different scaling param- 
eters a 7^ /? > 0, compact finite expressions involving 
Jacobi polynomials could be derived |T31 Eq. (4.6)]. 
In addition, there are infinite series representations 
expressing an overlap integral with different scaling 
parameters as an infinite series of overlap integrals 
with equal scaling parameters [131 Eq. (4.9)]. 

Explicit expressions of similar complexity could 
also be derived for nuclear attraction integrals ^13, 
Eqs. (6.4) and (6.5)] and in particular also for 
Coulomb integrals, which are six-dimensional inte- 
grals and which describe the electrostatic interaction 
of one-center charge densities represented by B func- 
tions T31, Eqs. (7.5), (7.7), (7.10), (7.13), (7.15), and 
(7.16)] 

Since Slater- type functions l with integral 
principal quantum numbers G N can be expressed 
by finite linear combinations of B functions [37[ Eq. 
(6.3)], the remarkably compact formulas for overlap, 
nuclear attraction, and Coulomb integrals of B func- 
tions can be used to write down analogous expres- 
sions for the corresponding integrals of Slater-type 
functions in a straightforward way. 

The extremely compact expression (|5.1|) was the 
basis for the construction of certain one-range addi- 
tion theorems by Filter and Steinborn ^39]. These 
addition theorems are expansions in terms of a com- 
plete and orthonormal function set based on the gen- 
eralized Laguerre polynomials. They converge in the 



mean with respect to the norm of the Hilbert space 
L^(M^) of square integrable functions. 

For the evaluation of the expansion coefficients of 
these addition theorems, which are essentially over- 
lap integrals, Filter and Steinborn only had to com- 
bine (j5.ip with the following expansion of general- 
ized Laguerre polynomials in terms of reduced Bessel 
functions [HI Eq. (3.3-35)] (see also Eq. (3.17) 
and Ref [23] on p. 2736]): 

e-^L(,")(2z) = (2n + a + l) 



E 



{-2YT{n + a + v+\) 
v\{n-v)\V{a + 2v + 2) 



k+i/2{z). (5.2) 
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The Fourier Transform of B 
Functions 



The remarkably compact explicit expressions for 
overlap, nuclear attraction, and Coulomb integrals 
published in [T31 [37] undeniably constituted a major 
achievement which provided considerable evidence 
that B functions indeed play a special role in the 
theory of multicenter integrals of exponentially de- 
caying functions. 

But many open questions remained. Filter had 
to develop in his PhD thesis _36, Sections 7-9] some 
very sophisticated, but also very complicated math- 
ematical techniques in order to obtain the remark- 
ably simple expressions for the integrals mentioned 
above. It was by no means clear whether and how 
these techniques could also be used profitably in the 
case of more complicated multicenter integrals. 

There was also the problem that in the Steinborn 
group the Fourier transform method, which already 
at that time had been considered to be one of the 
most or even the most important techniques for the 
evaluation of multicenter integrals, had for a long 
time a very bad reputation and was wrongly believed 
to be not very useful. This unjustified prejudice was, 
however, shattered when I showed that B functions 
possess a Fourier transform of exceptional simplicity: 

B-,(a,p) = [27^)-^/' J e-'P'-i3™,(a,r)d3r 
- (2A)V2 . 



(6.1) 

This is the most consequential and also the most of- 
ten cited result of my PhD thesis [HI Eq. (7.1-6) on 
p. 160]. Later, was published in ^ Eq. (3.7)]. 
Independently and almost simultaneously, (|6.ip was 
also derived by Niukkanen [HI Eqs. (57) - (58)]. 
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The exceptionally simple Fourier transform (|6.ip 
gives B functions a unique position among ex- 
ponentially decaying functions. It explains why 
other exponentially decaying functions like Slater- 
type functions with integral principal quantum num- 
bers, bound state hydrogen eigenfunctions, and other 
functions based on generalized Laguerre polynomials 
can all be expressed in terms of finite linear combi- 
nations of B functions (details and further references 
can be found in [42l Section IV] or [29l Section 4]). 

As remarked above, the derivation of the remark- 
ably compact explicit expressions for overlap, nu- 
clear attraction, and Coulomb integrals of B func- 
tions published in p!3l[37] was fairly complicated and 
required considerable mathematical skills. If how- 
ever, the Fourier transform (|6.ip is used, then the 
derivation of these expressions is almost trivial [43| . 

The Fourier transform (|6.ip of B functions is 
also a convenient starting point for the evaluation 
of more complicated multicenter integrals of B func- 
tions by a combination of the Fourier transform 
method with numerical quadratures. This was first 
done in articles by Trivedi and Steinborn [44], Gro- 
tendorst and Steinborn [45l |46] , Homeier and Stein- 
born [III UHl [49] , and Steinborn and Homeier [50] , 
and in the PhD theses of Grotendorst [5T and Home- 
ier [S2]. 

The key problem of this approach is that eval- 
uation by numerical quadrature can become pro- 
hibitively difficult because of the highly oscillatory 
nature of the integrands. In recent years, Safouhi and 
coworkers greatly enhanced the efficiency of this ap- 
proach by combining quadrature rules with suitable 
convergence acceleration techniques (see for example 
[53], [54]. [55]. [56] and references therein). This may 
well be the currently most promising approach for 
the evaluation of complicated molecular multicenter 
integrals of exponentially decaying functions. 



7 Bessel Potential Spaces 

On the basis of our current level of understanding, we 
can safely claim that the exceptionally simple Fourier 
transform (|6.ip is the most important property of B 
functions. It fully explains the advantageous features 
of B functions in the context of multicenter integrals 
of exponentially decaying functions. Nevertheless, it 
is certainly legitimate to wonder why it took such 
a long time to discover both B functions and their 
Fourier transform (|6.ip . 

For example, Geller [57l [58l [59] had studied in 



his work on multicenter integrals of Slater-type func- 
tions certain types of radial integrals, which contain 
the radial parts of Fourier integral representations of 
B functions as special cases. In my opinion, it is hard 
to understand that neither Geller nor anybody else 
realized that these special cases deserve special at- 
tention because they represent functions with highly 
useful features. 

I found out later that the Fourier transform (|6.ip 
of B functions had at least partly been known be- 
fore in a completely different context. Aronszajn and 
Smith [60] showed already in 1961 in connection with 
their work on functional analytic properties of so- 
lutions of certain inhomogeneous partial differential 
equations that [1 + with | e M" and a > 

is the n-dimensional Fourier transform of what we 
now call a reduced Bessel function with argument 
\x\ with X e W\ 

Their starting point was the following n- 
dimensional generalized modified Hclmholtz equa- 
tion: 

[l-vY\ix) ^ six). (7.1) 

Here, u, g: R" — > C are functions, x G M", a > 0, 
and is the n-dimensional Laplacian. 

Let us assume that u and g permit n-dimensional 
Fourier transformation. Then, 

u(^) = (27r)-"/2 / e-'^"'u{x)d''x, (7.2) 



m = (27r)-"/2y' c-'€-5(a;)d"a;, (7.3) 

and the Fourier transform of [l — V^] is the func- 
tion [1 -f- 1^1^]"/^. Thus, the partial differential equa- 
tion (17.111 can be reformulated as follows: 



(7.4) 



[l + |||2]a/2 

Consequently, for a given inhomogeneity g{x) the un- 
known function u(x) can be expressed as an inverse 
Fourier integral: 

m 



l{x) = (27r)-"/2 / gi 



The inverse Fourier integral on the right-hand side 
of (j7.5p can be expressed as the convolution integral. 



uix) = (2^)-"/^! G„(a;-y).g(y)d"y, (7.6) 

where the so-called Bessel potential Ga{x) is defined 
as follows: 

G^{x) ^ ^^-^^d-^. (7.7) 



Ernst Joachim Weniger: The Strange History of B Functions 



8: Outlook 



9 



Aronszajn and Smith [60l Eq. (2, 10) on p. 414] 
showed that Ga{x) is essentially a reduced Bessel 
function (compare also Samko's book [BTJ pp. 184 - 
185] and references therein): 



|a;|("-")/^j^(„_„)/2(|a;|) 
2-i+"/2r(a/2) 



(7.8) 



Here, K(^n_a)/2{\x\) is a modified Bessel function of 
the second kind [131 p. 66] 

As already discussed in Section [1] it is my convic- 
tion that science is now highly fragmented and com- 
munication between researchers from different scien- 
tific (sub)disciplines does not function particularly 
well. Bessel potentials provide convincing evidence 
supporting my claim. For example, mathematicians 
working with Bessel potentials seem to be completely 
unaware of the closely related Bessel polynomials and 
reduced Bessel functions, and in the literature of the- 
oretical chemistry Bessel potentials have apparently 
been completely ignored so far (an exception is my 
Habilitation thesis [521 P- 138]). 

8 Outlook 

Because of their remarkable simplicity in the coordi- 
nate representation, Slater-type functions have obvi- 
ous advantages. However, their molecular multicen- 
ter integrals are notoriously difficult. 

In contrast, B function have a comparatively 
complicated structure. Therefore, a more detailed 
look at their mathematical properties is required to 
understand why B functions have some highly advan- 
tageous features in connection with molecular mul- 
ticenter integrals. This does not mean that their 
molecular multicenter integrals are necessarily sim- 
ple in absolute terms, but they can in most cases 
be evaluated significantly more easily than the cor- 
responding integrals of Slater-type functions. 

Thus, the history of B functions, which is essen- 
tially a history of the influx of mathematical knowl- 
edge into theoretical chemistry, is well suited to study 
a question of more profound importance, namely 
the interaction between mathematics and theoreti- 
cal chemistry. 

I think that I can claim with some justification 
that the historical development of B functions has 
not been a straight path, but at best a not com- 
pletely random walk. Moreover, the available infor- 
mation has not always been utilized effectively, which 
obviously slowed down progress. Thus, communica- 
tion problems - not only between mathematicians 



and theoretical chemists, but also among theoretical 
chemists - have more or less been the rule rather 
than the exception. 

Since theoretical chemistry is such a highly inter- 
disciplinary field, future progress will depend very 
much on our ability of minimizing the detrimental 
effects of a breakdown of interdisciplinary and/or 
intradisciplinary communication. The history of B 
functions shows that this may not be so easy. 

The topics discussed in this article indicate that 
progress in theoretical chemistry depends very much 
on the influx of previously unknown mathematical 
ideas into theoretical chemistry. It would, however, 
be wrong to assume that communication between 
mathematics and theoretical chemistry would nec- 
essarily be a one way street. Problems from theo- 
retical chemistry in general and from the theory of 
B functions in particular can provide valuable inspi- 
ration for mathematical research. My own research 
provides ample evidence that this is indeed possible. 

During the work for my PhD thesis [TS], series 
expansions for multicenter integrals played a major 
role. Since it (too) often happened that my series ex- 
pansions converged slowly (see for example [HI Ta- 
ble II]), it was a natural idea to speed up conver- 
gence with the help of (nonlinear) sequence trans- 
formations. To the best of my knowledge, this was 
first done in 1967 by Petersson and McKoy |63|. Un- 
fortunately, I knew at that time only linear series 
transformations as described in the classic, but now 
outdated book by Knopp [5^. Unfortunately, these 
linear transformations turned out to be ineffective. I 
was completely ignorant of the more powerful non- 
linear transformations, which often accomplish spec- 
tacular improvements of convergence. 

My ignorance only changed when I did postdoc- 
toral work at the Department of Applied Mathemat- 
ics of the University of Waterloo in Waterloo, On- 
tario, Canada, where I - inspired by Jifi' Ci'zek - 
applied Pade approximants and continued fractions 
for the summation of divergent power series. After 
my return to Regensburg, I tried to apply nonlinear 
transformations also to slowly convergent series ex- 
pansions for multicenter integrals. In some cases, 
remarkable improvements of convergence were ob- 
served. 

In order to understand better the power as well 
as the limitations of nonlinear sequence transforma- 
tions, I also worked on their theoretical properties. 
As a by-product, I was able to derive several new 
transformations. The majority of these transforma- 
tions was published in my long article 65 , where also 
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efficient algorithms for the computation of sequence 
transformations as weU as theoretical error estimates 
and convergence properties are discussed. 

Later, I applied sequence transformations suc- 
cessfully in such diverse fields as the evaluation of 
molecular multicenter integrals of exponentially de- 
caying functions, the evaluation of special functions 
and related objects, the summation of strongly di- 
vergent quantum mechanical perturbation expan- 
sions, the prediction of unknown perturbation se- 
ries coefficients, and the extrapolation of quantum 
chemical crystal orbital and cluster electronic struc- 
ture calculations for oligomers to their infinite chain 
limits of stereoregular gwasz-onedimensional organic 
polymers. More information on my work both on 
and with sequence transformations and exact ref- 
erences can by found in my recent publications 

Thus, it is probably justified to claim that numer- 
ical mathematics and scientific computing ultimately 
profited via cross-fertilization from the convergence 
problems which I encountered during my PhD thesis 
and also later. 

Something similar can be said about the work of 
Safouhi. As already mentioned in Section [6l he com- 
bines numerical quadrature schemes with sequence 
transformations for the evaluation of the highly os- 
cillatory integral representations for multicenter in- 
tegrals. Over the years, Safouhi and his cowork- 
ers have worked hard and studied numerous differ- 
ent sequences transformations in order to optimize 
their computational approach. Initially, Safouhi was 
only interested in the evaluation of molecular multi- 
center integrals, but recently he could demonstrate 
that his techniques work very well also in the case of 
extremely pathological oscillatory integrals that are 
predominantly of mathematical interest |56| . 

Let me summarize. It is my conviction that the- 
oretical chemists interested in methodological ques- 
tions should actively seek the collaboration with 
mathematicians. Both sides have good chances of 
profiting from such a collaboration. 



A Notation and Terminology 

For the set of positive and negative integers, I write 
Z = {0, ±1, ±2, . . .}, for the set of positive integers, 
I write N = {1,2,3,...}, and for the set of non- 
negative integers, I write No = {0,1,2,...}. The 
real and complex numbers are denoted by M and 
C, respectively. The set of three-dimensional vec- 



tors r = [x, y, z) with real components j/, z £ R is 
denoted by R'^. 

For the commonly occurring special functions of 
mathematical physics I use the notation of Magnus, 
Oberhettinger, and Soni [T4] unless explicitly stated 
otherwise. 

For the functions of angular momentum theory 
- essentially (surface) spherical harmonics Y^{6,(j)), 
regular solid harmonics y]^{r), and Gaunt coeffi- 
cients (^3m3|£2W2|^iTOi) - I use the same notations 
and conventions as in [iTl Appendices B and C]. 

Fourier transformation is used in its symmetrical 
form, i.e., a function / : R^ C and its Fourier 
transform / are connected by the integrals 



Hp) = (2^) 



-3/2 



f{r 



(A.l) 



/(r) = (2^)-3/2y' e- P/(p)d3p, (A.2) 
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